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Abstract. In 1993 A. Ori [l] presented spacetimes violating the chronology 
condition in order to answer the question whether a time machine construction 
has to violate the weak energy condition or not. Later, in 2005 [2j, he con- 
structed a class of time machine solutions with compact vacuum core. Both 
classes include an interesting global structure and it is possible to obtain closed 
timelike curves. Besides we focus on the geometric structure, in particular sym- 
metries and geodesies, if feasible, and visualize several aspects. 



1. Introduction 

The spacetimes analyzed here, known as Ori-type spacetimes, were foremost de- 
scribed in [T] and [5] . The first paper presents a time machine model in which closed 
timelike curves (CTCs) evolve in a bounded region of space from a well-behaved 
spacelike initial slice. This slice S, just as the entire spacetime, is asymptotically flat 
and topologically trivial. In addition, this model fulfills the weak energy condition 
on S and up until and beyond the time slice, displaying the causality violation. 

The second paper uses particular vacuum solutions to construct time machine mod- 
els where the causality violation occurs inside an empty torus, constituting the core 
of the time machine. In that case the matter field surrounding the (empty) torus 
satisfies the weak, dominant, and strong energy conditions. 

Here we will comprehensively discuss the geometric and CTC structure of the class 
of Ricci flat chronology violating spacetimes given in [5] . The analysis for the other 
class discussed in [T] will only be given as a supplement, since — although it is the 
predecessor in some sense — the details are much harder to reveal because of the 
more complicated metric tensor, so here the research is still going on. 

The structure of other particular Ori-type spacetimes — pseudo Schwarzschild and 
pseudo Kerr — are discussed in [3]. 

Before starting with our main investigations we describe the important underlying 
time machine structure. 

2. Preliminaries 

Particles move through spacetime along causal curves, null curves are the trajecto- 
ries of massless particles whereas particles with non- vanishing mass are described by 
timelike curves. Suppose now that we are given a spacetime M containing CTCs. 
The physical interpretation would be that a particle moving once around a CTC 
would encounter itself in its own past. A similar argument can be adduced in the 



case of null curves. This is the reason why spacetimes with closed causal curves 
(CCCs) have been qualified as time machine spacetimes. 

While it is easy to construct spacetimes containing CTCs or CCCs, it is more diffi- 
cult to find examples of spacetimes that are — to some extent — physically reasonable. 
In [5], A. Ori presents a list, setting constraints on what should be classified as a 
realistic time machine model. 

One essential notion in General Relativity is the Cauchy development or domain 
of dependence of a subset A of a spacetime M. It is connected with the problem 
of formulating the Einstein equations as a well posed initial value problem and the 
concept of global hyperbolicity. 

From a realistic point of view it is desirable to obtain CCCs as the unevitable conse- 
quence of a Cauchy development. However, the interior of the Cauchy development 
of an achronal subset A is globally hyperbolic and, therefore, excludes any CCCs. 
The most one can hope for is that there are points in the boundary of the Cauchy 
development which lie on CCCs. This is the central idea of a time machine structure 
(TM-structure) for a spacetime. 

We will consider M to be a four-dimensional manifold and a metric tensor g on M 
with signature ( — h ++), such that the Lorentzian manifold (M,g) together with a 
fixed time-orientation is a spacetime. Let V(M) denote the set of all points p <E M 
that lie on some closed causal curve in M. We shall refer to V(M) as the causality 
violating region or the time machine. If p £ V(M), we say that causality is violated 
at p. The causality condition is said to hold on a subset U C M if V(M) n U = 0. 

We say that (M, g) has TM-structure if the following three properties hold for M: 

(TM1) There exists an open subset U of M on which the causality condition holds, 
(TM2) there is a spacelike hypersurface ,5^ contained in U, such that 
(TM3) for an achronal, compact subset ^ C 5^ the future domain of dependence 
of c to contains points in its boundary at which causality is violated, i.e., 
D + ('W) n V(M) 0. In this case the causality violating region V(M) is said 
to be compactly constructed. 

This seems to be a reasonable condition for any spacetime that may be viewed 
as a time machine model. Although there definitely are many more additional 
conditions one could impose, we will focus on TM-structure here. It should be noted 
that models with TM-structure can satisfy the weak, strong and dominant energy 
condition (see [5]) for the matter content of the spacetime if properly constructed. 

Example: Consider the smooth manifold M '■= S 1 x K 3 with the global coordinate 
system (</>, x, y, z) on M, cf> being a circular coordinate. Define the metric tensor on 
M by 

(1) g := —dcj) <E) dcf) + dx <£) dx + dy (g) dy + dz (£) dz. 

Then (M,g) becomes a vacuum spacetime, the so-called (four-dimensional) Lorentz 
cylinder if we require the timelike vector field to be future pointing. We see 
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immediately that all ^-coordinate lines are CTCs. Hence, every point of M lies on 
a closed timelike curve. 

It can be shown that V(M) = for any globally hyperbolic spacetime M (cf. [4|, 
Chapter 14). On the other hand, in the case of the Lorentz cylinder, we have 
V(M) = M. Mainly because of physical reasons (or everyday experience) we do not 
want to consider the Lorentz cylinder as a spacetime with TM-structure. 

The idea behind a TM-structure is the following: One can think of D + (^f) as the 
set of points p £ M that are predictable from ^ in the sense that no inextendible 
causal curve through q G D + ( c io) can avoid . Because the set int D + { c io) does not 
contain points at which causality is violated (for it is a subset of D + (ff) on which 
the so-called strong causality condition holds), the most one can hope for is that 
there are points in the boundary of D + (^) contained in V(M). 

As it can be seen from the definition of the TM-structure, spacelike hypersurfaces 
are an important entity in the context of TM-structures. Therefore, we include the 
following criterion. 

Proposition 1. Let f: M — > M be a smooth function on a Lorentzian manifold M 
and a G R be in the image of f . Then the sets 

(2) H t := {p G f-\a): ((grad/) p , (grad/) p ) > 0} 
and 

(3) H s := {p G /-»: ((grad/) p , (grad/)„) < 0} 
are hypersurfaces in M and H t is timelike, H s spacelike. 

Proof: The set U := {p G M : ((grad f) p , (grad f) p ) > 0} C M is open, hence a 
Lorentzian manifold in the usual manner. Because of the defining property of U the 
one form df\u is nowhere vanishing. Therefore, H t is a hypersurface in U, hence in 
M. If p G H t and u G T p H t , we have 

(4) 3((grad/) p , u) - u(f) - u{f\ Ht ) = 
since / is constant on H t . This proves the decomposition 

(5) T p M = TpH t © M(grad /)p 

and because (grad/) p is spacelike this implies that g\r v H t is of index 1, i.e., H t 
is timelike. If we consider H Sl almost the same arguments apply, except that now 
(grad f) p is timelike. Hence T p H s is a spacelike subspace of T p M, i.e, H s is spacelikc. 

□ 

In coordinates x 1 , . . . , x n the differential df of / is df — -§£tdx z and 

(6) (grad /, grad f)=g^^L^l. 

Using this formula and applying the preceding proposition to the coordinate neigh- 
borhood, one obtains: 
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Proposition 2. If x ,. . . , x n are local coordinates of a Lorentzian manifold M , the 
coordinate slices x k = const for 1 < k < n are 

(i) spacelike hyper surf aces if g kk < 0, 

(ii) timelike hypersurfaces if g kk > 

for all points p in the coordinate neighborhood. 



3. A Class of Ricci Flat Chronology Violating Spacetimes 

Here we investigate spacetimes described in [2] and we start with the smooth man- 
ifold M := M 3 x S 1 and introduce coordinates (T, x, y, <j>) on M, where T, x and y 
are natural coordinates on IR 3 and <fi is a circular coordinate on 8 . For any smooth 
function /: R 2 x S 1 -> 1 we define the metric tensor g on M by 

(7) g := -dT ® s d(f> + dx ® dx + dy ® dy + (f(x, y, 0) - T)d<j) ® d(f>, 

with ® s being the symmetrized tensor product. Emphasizing the role of /, we 
write Mf for the semi-Riemannian manifold (M,g). The periodicity of / in its 
third argument guarantees that g is well defined on all of Mf. Since the matrix of 
the component functions of g has detfgy] = — 1, the metric tensor is non-degenerate 
and has Lorentzian signature. 

The Ricci tensor of Mf is given by 

1 fd 2 f d 2 f\ 

(8) Ric -- 2 (^ + ^) d ^^ 

so we assume that / satisfies the relation 



(9) fxx ~t~ fyy 

to get Ricci flat spacetimes Mf. As the timelike unit vector field 

(10) d^ + ±(f-T+l)d T 
shows, Mt is time orientable. 



3.1. General Properties. Before we make an explicit choice for / and embark on 
a more detailed study of a special example of the described class, we shall consider 
the general case. 

3.1.1. Closed Timelike Curves in Mf. Of course, it is the topological factor S 1 in 
Mf and the metric component g^ = f — T that is responsible for the appearance 
of CTCs. Given T, x and y the curves 

(11) [-7T,7T] -> M, S I— > "fT,x,y (s) •= (T, x,y, s) 

are closed and timelike provided 

(12) f(x,y,s)-T<Q for all sG[-7r,7r]. 
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Depending on the explicit form of /, this is a condition on the coordinates (T, x, y) 
determining a non-empty subset of IR 3 that consists of points sitting on CTCs of 
the form ( 11 ). By Prop. [2] the hypersurfaces Ht defined by T = const are 



, timelike] f > o, 

(13) rl if 9 = T — f is I 
spacelike J < 0, 

respectively. Therefore, any timelike curve in Ht, in particular CTCs, must be 
contained in that part of Ht where T — f > 0. The same result may also be 
obtained by calculating 5(7', 7') directly with the condition dT("f') = 0. We will 
see later for one choice of / how these two ingredients — the region of CTCs and the 
causal character of Ht — allow for a TM-structure of Mf. 

3.1.2. Killing Vector Fields of Mf. Let C be the Lie derivative, such that the Killing 
equation for a vector field X is given by Cxg — 0. Furthermore we denote the 
algebra of Killing vector fields on some open subset U C M of a manifold M by 
i(U). We will make repeated use of the following proposition, the proof of which 
can be found in [TU], p. 61. Here R denotes the curvature tensor of M. 

Proposition 3. If X 6 i{M) is a Killing field, then 

L x D {k) R = for all k € N, 

i.e., the Lie derivative in the direction of X of any covariant derivative of the 
Riemann curvature tensor vanishes. Let L^ denote the linear system (C X D^R) P 
for some p G M and k £ N. Then 

dimi(Af) < dimkerLfc. 

For arbitrary /, depending on all three coordinates x, y and <fi, there is no reason 
to expect an abundance of non-trivial, linearly independent Killing vector fields on 
Mf. There is only one obvious and non-tivial Killing vector field, which is actually 
not even globally defined. 

Proposition 4. For arbitrary f there is only one local Killing vector field K on 
Mf. In coordinates, 

(14) K = e~^d T - 

Proof: First of all, note that K cannot be extended to all of Mf. Formally not 
entirely correct, we may say that this is due to the fact that the component function 
e~i of K is not properly periodic in <p. If 4> takes values in (— 7r,7r) the problem 
arises at points p € Mf which are not covered by the coordinate system (T, cc, y, <p). 
In order to have an atlas of Mf at our disposal, we use a second coordinate system 
(T, x, y, ip) with tjj mapping into (0, 2ir). We shall refer to the first coordinate system 
as A and to the second as B. The transformation formula 

[cf) if <(j)< n, 

( 15 ) = \ J. o 

\ 6 + 2tt if -7T < 6< 0, 



implies the coordinate expression 

{e~ 2 d T if O<0<tt, 

* 
e-i^dr if -7r<</><0, 

for K in system B. li K were extendible to all of Mf, using B, the two limits 

(17) ]imK p = e-*dr 

and 

(18) lim K p = lim e" _7r <9T = e~~dr 

would necessarily be the same. Next we verify that K is indeed a Killing vector 
field. Generally, for a vector field X :— Gdr + Hdx + Pdy + Qd,/, the Killing equation 
Cxg — is equivalent to the following system of PDE's: 

Qt = 
H t -Q x = 
Pt -Q v = 

Py=0 

H x = 

P X + Hy = 

G T - fQr + TQ T + Q = 
G x ~ fQx + TQ X — — 

Gy — fQy + TQy — Py = 

-2G + 2Q / + QU ~ 2TQ^ + Hf x + Pf y - G = 

The assumption H = P = Q = and the supposition that G depends only on <fr 
reduces this system of PDEs to the simple ordinary differential equation 

(19) G + 2G = O. 

Hence G(<p) = and if is a Killing vector field. It remains to show that there 
is a choice for /, such that there is no other solution to the Killing equation which 
is not a constant multiple of K. This is the point where Prop. [3] comes in. If we 
choose 

(20) f(x,y,(f>):= encash 

the equation CxDR = evaluated at (T,x,y,(j>) = (1,1,1,0) becomes a linear 
system with kernel of dimension 1 . □ 

The fact that K cannot be extended to a global Killing field is due to the topology 
of Mf. If we consider the universal covering k: IR 4 — > Mf of Mf, assign to ]R 4 the 
pull-back metric k*g and use global coordinates (T, x, y, z) on K 4 , the vector field 

(21) K:=e-id z 

is a global Killing field on M 4 , irrespective of the explicit form of /. Changing the 
^-coordinate to the periodic coordinate 4> implies loosing the global Killing field K. 



3.2. One Special Choice for /. We are now going to consider a special example 
of the class of spacetimes described so far by choosing / to be 

(22) f(x,vA)--=^ 2 -y 2 ), 

where a is a positive constant. For notational convenience we drop the subscript / 
in Mf from now on. In order to obtain a coordinate system which is better suited 
to the description of the region of CTCs in M we use the transformation: 

(23) M 3 x S 1 -> R 3 x S\ (T,x,y,(j>) >-> (t,x,y,</>) 
with 

(24) t:=T _^_ y 2 ) + ep 2 

Here, e is another positive constant and p 2 := x 2 + y 2 . It is clear that (t, x, y, <fi) is 
a coordinate system. A short calculation yields 

(25) g = —dt ® s d(f> + (2e — a)xdx ® s dcf) + (2e + a)ydy ® s d<j> 

+ dx ® dx + dy ig> dy + (ep 2 — t)d(j) ® d</> 

in the new coordinates. 



[ spacelike 


for 


t < ep 2 , 


null 


for 


t = ep 2 , 


[ timelike 


for 


t > ep 2 . 


are spacelike. 


At t = 


the curve 70,0,0 is a closed null 



3.2.1. CTCs in M. For fixed values of t, x and y the curves Jt,x,y(s) = (t,x,y, s) 
are closed and 



(26) 



geodesic, as will follow from later results. For a given value of t all curves jt,x,y 
lie in the hypersurface H t given by t = const. The causal character of H t can be 
inferred from 

(27) g u =t+ (2e - a) 2 x 2 + (2e + afy 2 - ep 2 . 
If we suppose 

(28) (2e + a) 2 < e, 

then 7f t is spacelike throughout for t < 0. At < = 0, H t is spacelike except at the 
closed null geodesic x = y — 0. In the region i > 0, 

{spacelike if (x,y) lies outside E t , 

timelike if (x,y) lies inside E t , 

where E t is the ellipse in the x-y plane described by g u = or 

(30) * = [e - (2e - a) 2 ]a; 2 + [e - (2e + a) 2 ]y 2 . 

Thus any CTC in iJ t must be entirely contained in E t (see Fig. [T]). 

In [2], Ori describes in some detail how to construct a spacelike hypersurface 5? in 
the region {p e M : < 0} of M that contains a compact subset If, such that 
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Figure 1. A section <j> = const in H t and the circle of radius t 
in the x-y plane filled with CTCs of the form Jtx,y ( we have set 
e = 1). The growth of this circle with t is indicated by the arrows. 
The ellipse E t confines the region of CTCs in H t . 



the closed null curve 70,0.0 is contained in £H 
spacetime with TM-structure. 



making M a candidate for a 



3.2.2. Killing Vector Fields of M. We have already seen that due to the topology 
of M the Killing field K, which exists on M independently of the explicit form of 
/, is not a global Killing field. We will be accompanied by this problem throughout 
this section and, therefore, we concentrate our efforts on the part U of M covered 
by the coordinates (i, 2, y, </>), where <fi € (— 7r, tt). However, since the expression for 
/ does not contain the coordinate <j> we immediately conclude that the vector field 
dtf, is a global Killing field on M; it will be the only one. 

Proposition 5. The dimension of i(U) satisfies 



(31) 



dimi(Z7) < 6 



and Kq — and K^ = e 2 dt are Killing fields on U . 



Proof: We use Prop. [3] to get the upper bound for dimi(U). At the point 
p = (1,1, 1,0) the linear system CxR = consists of the following four linearly 
independent equations: 



{L X R)4242 = -^X 4 - (2c - a)aX U2 + aX h4 

a a 
{L x R)i3i2= 2 Xl '' 2 ~ 2 ~~ a ) Xl ;3 - aX 2-3 

Hence, the result follows. Since dt — &r, the vector field K$ is a Killing field because 
of Prop. |4j □ 



The indices of K§ and K% in the foregoing proposition indicate that indeed 



dimi(P) — 6. In order to find the other four Killing fields we have to write down 
explicitly the equations given by the Killing equation Cx9 = 0. If we set 

(32) X := Gd t + Hd x + Pd y + Qd^ 
these are: 

(33) Q t = 

(34) H t + (2e - a)xQ t - Q x = 

(35) P t + (2e + a)yQ t - Q y = 

(36) -G t + (2e - a)xH t + (2e + a)yP t + {ep 2 - t)Q t - Q = 

(37) H x + (2e - a)xQ x = 

(38) P x + (2e + a)yQ x + H y + (2e - a)xQ x = 
-G^ + (2e - a)i? + (2e - a)x.ff x + (2e + a)^ 

(39) +(e,o 2 - i)^ + + (2e - a)a;Q0 = 

(40) P y + (2e + a)yQ„ = 
-G y + (2e - a)xffy + (2e + a)P + (2e + a)yP y 

(41) +(e j o 2 -i)Q y +P^ + (2e + a)yQ = O 

-G - 2tQ + 2e(a:iJ + j/P) + 2ep 2 Q 4> 

(42) +2(2e + a)yP + 2(2e - a)xiJ - 2G = 



The following two assumptions will simplify this system: 



• Q = 0; 

• H and P are functions of (f> only. 
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Then eqs. (33l-(35l, (37 1, (38 1 and eq. (40 1 are trivially satisfied and we are left 
with: 

(43) G t = 

(44) -G x + (2e - a)H + P = 

(45) -G y + (2e + a)P + P = 

(46) -G + 2e{xH + yP) + 2(2e + a)yP + 2(2e - a)xH^ - 2G = 
At this point there are two similar cases. 



1. Suppose P = 0. The last three equations simplify to: 

(47) G„ = 

(48) —G x + (2e - a)# + P = 

(49) -G + 2exH + 2{2e-a)xH (j ,-2G 4> = 



Because of (43 1 and (47 1, G does not depend on t or y. Since P is a function of . 



only we deduce from ( |48| that 

(50) G(x,^) =x(F^) + (2e-a)F(<£)), 

where we have set a possible integration constant equal to zero. Substituting this 



into (|49j results in the ODE 

(51) 2H^ + H <i ,-aH = 0. 
The ansatz H{<f>) = e" x ^ leads to the polynomial equation 

(52) 2A 2 - A - a = 
with zeros 
(53) 

and the solutions to ( |51| are 

(54) :=e- c " t ' and H 2 {4>) := e'^ . 
For G we calculate 

(55) Gi(x» := (2e - a - o)a;e~ Q0 and G 2 (x, 0) := (2e - a - /3)xe _w . 

With these choices for P and G all 10 components of the Killing equation are sat- 
isfied. 



a := \ (l + Vl + 8a) and /3 := - (l - Vl + 8a) , 



H. Suppose H = 0. From (44 1, (451 and (46) we now get: 

(56) G x = 

(57) -G y + (2e + a)P + P = 

(58) -G + 2eyP + 2(2e + a)yP - 2G = 
Reasoning as in the first case leads to 

(59) G(y,cf>) = y(P^ ( f>) + (2e + a)P( ( f>)) 
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and an ODE for P: 

(60) 2P 00 + P + aP = 
The solutions for this equation depend on the value of a. 

A. For < a < | they take the following form: 

(61) P n (<t>) ■■= e-^ and P 12 (0) := e~ v + , 
where 

(62) /i := ^ (1 + VI - 8a) and i/ := 1 (l - Vl - 8a) . 
The corresponding expressions for G are 

(63) G 11 (y,cl>):=(2e + a-fi)ye-^ and G 12 (y, 0) := (2e + a - 
P. If a = | , the solutions are 

(64) P 2 i(0):=e-i; P 22 (<f>) := <j>e~ % ; 

(65) G2i(l/,0):=(2e-^)i/e-i; G 22 (y, 0) := ((2e - ^ + l)ye" 
and finally: 

C. For a > | we have 

(66) P 3 i(0) := e-i sin A(0), P 32 (0) := e~* cos A(0) 
and 

(67) G 3 i(y,4>) := Q\/8^TcosA(0) + (2e + a- ^sinA^ye-s 

(68) G 32 (y» := ^- * V8a - 1 sin A(</>) + (2e + a - J) cos A(<f>)j ye - * 
with A(</>) := |V8a -10. 

We summarize: 

Proposition 6. A basis for i(U) is given by the following Killing fields: 

K^e-^i^e-a-^xdt + dx] 
K 2 = e-^[{2e-a-fi)xd t + d x } 
K 3 = G tl d t + Pad y 
K A = G l2 d t + P i2 d y 

K e = dj, 

li 



Here, the Go? and are the functions defined above and 




Note that, as in the case of the Killing field K in Prop. [4] none of the Killing fields 
K\ - K§ can be extended to all of M. 

3.2.3. Geodesies of M. As we shall see, with the given special form of /, it is 
possible to solve the geodesic equations on M analytically. We express the geodesic 
7: I — > M defined on some interval I around zero as 7(5) = (t(s), x(s), y(s), 0(s)), 
where s is a proper time parameter in the case of a timelike geodesic 7 and an affine 
parameter if 7 is lightlike or spacelike. Then the geodesic equations read: 

(69) 0-^0 2 = O 

(70) x-^ax^^O 

(71) V+\ayi 2 = Q 

i+i<fi- (2e - a)x 2 - (2e + a)y 2 - 2e(xx + yy)(j> 

(72) +- (t + (a 2 - 2ea - e)x 2 + (a 2 + 2ea - e)y 2 ) 2 = 
Furthermore, we have the condition 

(73) - 2i4> + x 2 + y 2 + 2{2e - a)xx<j) + 2(2e + a)yy4> + (e(x 2 + y 2 ) - t)<p 2 = k, 
where k = —1,0, 1 for timelike, lightlike or spacelike geodesies, respectively. 



The solving process consists of two major steps. In the first step we solve eqs. ( 69 ) 



( 71 1 , which are independent of t, and in the second step we use ( 73 ) to determine 



an expression for t(s). 



Step 1. Use the Killing field K§ or directly integrate (69 1 to get 
(74) e~H = A 



for an arbitrary constant A. Note that (74 1 is equivalent to (69 1. The general 



solution for ( 74 ) in the case of A 7^ is 



(75) fts) = -21n (-^{s + B)Y 

where B is another constant. The case A = will be treated seperately later on. 
Since we would like 7 to be defined on some interval around zero, we assume B =/= 0. 
Because A has the same sign as cb(0), we see that <j) is on ly defined on 

(i) (-oo, -B) if 0(0) > 0, or 

(ii) (-5,00) if 0(0) < 0. 



In either case 7 is incomplete. By traversing 7 in the opposite direction, i.e., by 
setting -y(s) := j(—s), it is always possible to obtain the case (f>(0) < 0. Then 
s + B > on / and since this will simplify the calculations to come, we make the 
assumption 



(76) 



(f)(0) < 



B > 



A < 0. 



If we substitute the result ( 75 1 into ( 70 ) , we arrive at 

(77) *-(^By x = - 

The solution of this equation is 

(78) x(s) = C(s + B) 2a + D(s + Bf 13 



with a and (3 defined in ( |53| and two further constants C and D. Similarly, (71 1 
becomes 

2a 



(79) 



y- 



ry = o, 



(s + B) 2 ' 

but the solution to this equation depends on the value of a. With constants E, F 
and the function 

1 



k(s) := ^V8a- lln(s + 5) 

we have (for the constants \x and ^ confer (|62|): 
1 



(80) 
we h 

(81) 1! o 



y( S ) =E(s + B) 2fl + F(s + B) 2u 



(82) 
(83) 



1 



1 



y[s) =Vs + B(E + F\n(s + B)) 



a > - : y(s) =Vs + B [Esin(/t(s)) + Fcos(k(s))] 
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Step 2. To find an expression for t(s) is more complicated. First, we rewrite the 



metric condition ( 73 ) as 

(84) 

where 



k + 2i<f> + t<j) 2 = h, 



(85) h := i 2 + y 2 + 2(2e - a)xx<j) + 2(2e + a)yy<t> + (e(x 2 + y 2 ) - t)<j) 2 . 
This equation is easily integrated. One obtains 



(86) 



t(s) 



1 



(fc - ft(s)) ds + G 



(s + B), 



G being another constant. For each of the three possibilities for a we calculate h(s) 



and perform the integration in ( 86 1 . Since the explicit steps of the calculation are 



rather long but straightforward we suppress them and simply state the results. 



A. In the case of < a < | we have 



(87) 
with 



t(s) = i(s) + u 3 (s + B)^ +u 4 (s + Bf v 



t( s ) := -s(s + B)+ G(s + B) + m(s + Bf a + u 2 {s + B) 



±0 
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and the abbreviations 



(89) 

B. Here a 
(90) 
where 
(91) 



C 2 D 2 
ui := — (2e - a- a), u 2 := — (2e - a - 0), 

E 2 F 2 
u 3 := — (2e + a - fi), u 4 := — (2e + a - i/). 

g and 

t(a) = F(s) + («! + w 2 ln(s + B))(s + B) ln(s + B), 



"1 



lF(B(8e-i)-F) and v 2 := ±F 2 (8e - \). 



Finally, there is: 



C. For a > | the solution takes the form 

(92) i(s) = i(s) + [w 2 cos 2 (k(s)) — W\ cosk(s) sin«;(s)] (s + B), 

where the constants w\ and u>2 are 

V8a - l(E 2 - F 2 ) - (4(2e + a) - 1)BF , 



(93) 



w 2 



- - 2(2e + a) ) (£ 2 - F 2 ) - V8a - IEF 



In order to complete our programme we need 



Proposition 7. Any curve 7 satisfying the first three geodesic equations (69 I - (711 



and the metric condition ( 73 I is a geodesic 



Proof: We have to check that the geodesic equation for t is satisfied by 7. If we 
set 



(94) 



7" = a^d t + a^d x + a^d y + o^^, 



then equations (69l-(71l give a 



Jx) _ „( B ) 



condition (7', 7') = A: w.r.t. s results in 
(95) (7",7') = ^ 



l (4>) — 0. Derivation of the metric 



o<*> = 0. 



As we know from (74 1, the function is either always zero or never zero. Thus, if 



4> is never zero, then — and 7 is a geodesic. In a moment we will see that the 
case </> = will lead to a geodesic as well. □ 



The case A = 0. This condition is equivalent to <f> — and therefore we immediately 
conclude 

(96) x(s) = c\s + Co and y(s) — d\s + do- 

The metric condition now reads 



(97) 



cf + df = k. 
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Hence for a causal geodesic, we have k = C\ = d\ = 0. Then the geodesic equation 
for t yields t — 0, thus 7 is of the form 

(98) j(s) = (ks + l ,x ,y a ,<p ), 

for constants and ii, i.e., linear parametrizations of ^-coordinate lines are null 
geodesies. If k = 1, we find that 

(99) 7(s) = (/ 2 s 2 + ^is + /o,cis + c o ,dis + rfo,0o), 
where 

(100) l 2 :={2e-a)c\ + {2e + a)d\ and + = 1, 
are spacelike geodesies. 

3.2.4. Discussion. To simplify later expressions, we introduce the notation 

-00, if A > 0, 

(10J) isgn(A) := { if A = 0, 

00 if A < 0, 

for any real number A. As we see from the solutions, the constant B determines 
the domain of 7. Since we assume <p(0) < 0, we always have 7: (— B,oo) — > M. 
The causal character of 7 does not have any influence on the behavior of the three 
component functions tfi, x and y. Hence the following holds for any geodesic. 

Behavior of <fi. 

From the result ( 75 1 obtained above we determine 



(102) <Hs) = 2 —, B = -J-, i = <^(0)e-^. 

V ' W s + B 1 0(0)' VV ' 

Hence on I — (— B, 00) the function ip is strictly decreasing and its asymptotic 
behavior is (cf. Fig. [2]) 

(103) lim (j)(s) = +00 and lim 4>(s) — —00. 

s\ — B s— ¥00 

However, we must not forget that (f> is a circular coordinate, i.e., geometrically im- 
portant is the behavior of s i-> g(<f>(s)) := (cos 4>(s), sin 4>(s)) € S 1 . In particular, 
this means that as s approaches ~B from above, 7 circles infinitely often around 
S 1 in M = R 3 x S 1 . 



Behavior of x. 



Since a > and j3 < 0, we conclude from (78 1 that 



(104) lim x(s) — isgn(Z?) and lim x(s) = isgn(C). 
The initial values a;(0) and i(0) determine the parameters C and D by 

C = 1 - - ( x(0) - ^x(0) 

(105) VTTTaB^ { B 

1 / . 2a 



\/TT8^B 2 / 3 
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_ ^(o) - ^(O^j 



Figure 2. Some examples of the geodesic component functions <f> 
and x. On the left the parameter values are A = — 1, B = 1 (blue), 
A = -2, B = 1 (red) and A = — 1, B = 2 (green). For £ = -1 
the asymptote is shown. With a = ^ (right), the function x is 
displayed for C = D = 1 (blue) and C = — 1, D = 1 (red). 

Behavior of y. 

Here again, we have to distinguish the already well known cases for a. 



A. Because of < 4z^ < 1 < 4/i, (81 1 implies 

lim y(s) 



and 



isgn(£) 
isgn(F) 

The parameters E and F are related to the initial values y(0) and y(0) by formulas 



(106) lim y(s) = 

s~\—B 



if £^0, 
if E = 0. 



completely analogous to ( 105 1 



5. It follows from (82 1 that 



(107) lim y(s) = 

s\—B 



and lim y(s) 

s— > 00 

The relevant parameter f can be obtained by 



isgn(F) 
isgn( J B) 



(108) 



F 



B j/(0) 



y(o) 

2B 



C. In this case we refer to ( 83 1 . Still 
(109) 



if F^0, 
if F = 0. 



lim w(s) = 0, 

S\-D 

but for s — > 00, the function y(s) does not converge, in general. Instead it oscillates 
consecutively between the values F, E, —F and — E multiplied with y/s + B. 



Behavior of t. 

It is the form of t that shows the causal character of geodesies. Since the expression 
for t is different in each of the three cases induced by the parameter a and does 
generally include all parameters from B to F, there are too many possibilities of 
asymptotic behavior to plot examples for all of them. Therefore, we are going to 

16 



Figure 3. The geodesic component function y. For a = ^ (left) 
we chose E = F = 1 (blue) and E = — 1, F = 1 (red). On the right 
a = | and the parameters are E = F = 1 (blue) as well as E = 1, 
F = -1 (red). 

6- 



y(s) 0- 



6 1 2 ' 3 4 ' 5 

S 

Figure 4. When a = 1 we get oscillating behavior for y(s). E = 1. 
F = 2 is in blue. E = —1, F = 2 in red and £ = 2, F = — 1 green. 




illustrate the behavior of t by some examples focussing on causal geodesies and 
mention some properties as we go along. 



A. In this case the expression for t is given by (87 1. The inequalities 
(110) 4/3 < < Au < 1 < Ajj, < 2 < 4a 

imply the limit 



(111) 



lim t(s) — isgn(u2) oo. 

s\,—B 



Using the same inequalities we also get 



(112) 





is? 


m(ui) 


if 


til + o, 




isf 


rn(fc) 


if 


mi = 0, k i= 0, 


lim t(s) = < 

s— >oo 


is{ 


m(u 3 ) 


if 


u\ = k = 0, W3 7^ 0, 




isj 


rn(G) 


if 


ui = 1*3 = fc = 0, G ^ 0, 






m(w 4 ) 


if 


iti = U3 = k = G = 0. 



B. Analyzing (J90|) , we come to the conclusion that the asymptotic behavior of t for 
s 



-B is the same as in case A. Furthermore, 



(113) 



lim t(s) 



isgn(iti) if itiT^O, 

isgn(fc) if iti = 0, k ^ 0, 

isgn(w 2 ) if ui = k = 0, v 2 ^ 0, 

isgn(wi) if u i — k — v-2 = 0, v\ ^ 0, 

isgn(G) if ui = v 2 = v 1 = k = 0. 



In the following figure, where we illustrate cases A and B, we have u\ > and m 2 > 0. 




S s 

Figure 5. The function t for a = ^ (left) and a = | (right). In 
either case the red curve is a lightlike geodesic, whereas the blue 
one is timelike. All parameters, except a, are equal to 1 in both 
plots. 



C. Here we have to deal with (92). The behavior of t as s tends to —B from above 
is the same as in the foregoing two cases. This is also true for lim^oo t(s) except 
if Mi = k = 0. In that case, t does in general not converge as s tends to infinity. 

Thus in Fig. [6] the curve in red is the i-component of a lightlike geodesic with 
Mi = 0. For the blue ^-component of the timelike geodesic shown we have u\ = 
as well, but since here k = —1, there is no oscillation and lim^oo t(s) = — oo. 




10 20 30 40 50 60 70 



Figure 6. The function t for a lightlike geodesic (red) and a time- 
like geodesic (blue) in the case of a > |. Here, all parameters, 
including a, are equal to 1. 

4. Ori's Spacetime of 1993 

Our exposition here is mainly based on |TJ , but we also refer to [7] and [5] for further 
details about the physical relevance of this spacetime. 



4.1. The Manifold. Topologically the manifold has the form M := M 4 . In order to 
specify the metric we start with standard coordinates (t, x, y, z) on M and introduce 
polar coordiantes on the planes of constant t and z according to the usual formulas 



(114) 



y 



r cos i 



r sin < 



with r € (0, oo), <f> <E (— it, it). 



Interpreting as a circular coordinate, we obtain the alternative coordinate system 
(t, r, <f>, z) on M \ S, where S := {p £ M : x — y — 0}. Away from S we will use 
these coordinates from now on unless otherwise mentioned. The metric then reads 
(115) 

" -1 ahrt 

1 -bhr(r-r Q ) 



\9 



ahrt 



-bhr(r-r ) r 2 (l + h 2 (b 2 p 2 - a 2 t 2 )) -bhrzO -bhrz 1 



Here a, b, r > are constants, p(r, z) 2 := (r — r Q ) 2 + z 2 and h 



is a function 



of class at least C 2 depending on p satisifying the following conditions: 

(i) < h{p) < 1 for all p, 

(ii) h(p) = 1 for p < 0, 

(iii) h(p) =0 for p > d, 

(iv) h'(p) < for < p < d. 

For the fourth constant d we stipulate < d < ro. The sets {p £ M: t(p) = 
to,p(p) = pq < tq} are tori (in particular this is the case for p = d ) and when 
talking about such tori as being of the form p = const we tacitly assume t = const. 
The function h confers a special role to the parameter d because by means of h the 
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Figure 7. Illustration of the torus T. Outside T the metric is 
equal to the Minkowski metric. 

metric g reduces to the Minkowski metric 



outside the torus T := {p e M: p(p) = d}, which at the same time shows that g 
is well defined on M since the gap S is filled in by setting g to be the Minkowski 
metric on S, too. The component functions of g are at least C 2 everywhere, and 
we calculate det^j] = — r 2 . Hence, g is non degenerate. One can easily verify that 
the vector fields 

£?i := <9 t , E 2 :=d r , E 3 := ahtd t + bh(r - r )d r + + bhzd z , E A :=d z 

define an orthonormal frame on the range of the coordinates (i, r, <fi, z), and the rela- 
tions (Ei, Ei) = — 1 and (E k ,E k ) = 1 for k ^ 1 show that g has correct Lorentzian 
signature (throughout M). Everything presented here is independent of the explicit 
form of h, provided the conditions (i)-(iv) listed above are satisfied. However, there 
are circumstances when more information about h is needed. Calculating numeri- 
cal values of certain quantities or the explicit form of geodesies are such examples. 
Algebraically easy to handle is Ori's choice for h which is 



In this case one encounters the drawback of h being continuously differentiable up 
to order 2 only. A different choice which results in a smooth version would be 



(116) 



— dt® dt + dr ® dr + r 2 d(j) <E> d<j> + dz ® dz 



(117) 




(118) 



1 



P<0, 

< p < d, 

p>d, 




Hp) = { 




where 



t > 0, 



t < 0. 



Here h is a typical cutoff function as used in differential geometry. 
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4.2. CTCs in M. For any given real number t let H to denote the hypersurface 
t _1 (to) in M. We are interested in characteristics of the set 

Mcv(t) ■= {p € Ht : 37: I —> H t , p G 7(1), 7 closed and timelike}, 

where the subscript is meant to stand for chronology violation. First consider the 
curves 7: [— 7r,7r] — > H t , s ^ (t,r,s,z) for constant values of t,r and z satisfying 
p(r, z) < d (which implies h(p) 7^ 0, a fact we will need later). We have 

(119) (7', 7') = (d*, 84) =r 2 (l + h 2 (b 2 p 2 ~ a 2 t 2 )) 

which is independent of the curve parameter s. To check for the sign of this expres- 
sion we need to look at 

(120) l + h 2 (b 2 p 2 -a 2 t 2 ) = a 2 t 2 = b 2 p 2 + hT 2 . 

The function [0,d) — > [1, 00), p M> b 2 p 2 + h~ 2 is strictly increasing (in view of the 
properties of h) and surjective (hence bijective). Therefore, ( |120| has a unique 
solution p t iff \t\ > i. The function [^,00) — > [0,d),t i-> p t is likewise strictly 
increasing. Given these results we deduce from ( |119 1 and (120) that for at > 1 
(now we index 7 by two "parameters" r, z) 



timelike if p(r,z) < p t , 

(121) 7 r , z is <(null if p(r,z)=p u 



f 

^ spacelike if p(r, z) > . 

Throughout the range |t| < ^ the quantity ( 119\ above is always positive and 7 is 
spacelike, irrespective of the values of r, z. In conclusion we know that 

{p £ H t : p(r,z) < p t } <Z M C v{t). 

But we are also able to confine Mcv(t) as follows. Calculation gives 

(122) g tt = h 2 a 2 t 2 - 1 

and by Prop. [2J H t is a spacelike hypersurface for \t\ < ~. The situation changes 
when t > -. Then we have 



(123) H t is 



{timelike for p < p t , 

spacelike for p > p t , 



where p t is determined by h(p t ) = ^j, i.e., /J t is the (unique) zero of (122). Com- 
paring the equations for p t and p t yields the first part of the following inequality 
and the second part results from the properties of h: 

(124) p t < p t < d. 

In any given hypersurface H t the region p > p t is thus free of timelike curves which 
implies Mcvif) C {p G M: p(p) < pt}. In particular, all CTCs in H t , including 
the ones in the torus p = p tl are contained in the torus p = p t . This situation is 
illustrated in Figure [8] 

Note that for all t, viewed as subsets of H t , the set Mcv(t) is contained in the 
compact set T. 
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Figure 8. For fixed t > - this is a section of the torus T containing 
the x and z-axis (the t and y-coordinates are suppressed). The 
growth of the two smaller tori p — p t and p — p t with t is indicated 
by the outward pointing arrows. 

4.3. Tipping of Lightcones. The fact that the curves / y rz from above are CTCs 
filling a torus of growing radius p t , implies that the tangent vector 7^. z = becomes 
timclikc at points inside the torus T, where it was null or spacelike before (i.e., for 
smaller values of t). This is reflected by the following argument, where we place 
ourselves at r = ro, z = and consider a tangent vector of the form X = ad$ + f3d z 
with real coefficients a, f3. We have 

(125) (X,X) = p 2 + a 2 rl{l- a 2 t 2 ), 

which is positive for small values of t, zero at a certain to and negative for t large 
enough. Correspondingly the causal character of X is spacelike first, lightlike at 
to and then timelike. Using the concept of lightcones, X is initially (i.e., t small) 
situated outside the lightcone at the corresponding point of M, at to it lies on the 
lightcone and can finally be found inside the lightcone when t has grown sufficiently. 
This process of tipping of lightcones shows geometrically how CTCs emerge and is 
depicted in figure [9] (from [3]), where we can see the set {(t,r,(/),z) G M : z = 0}. 
The coordinate t increases in the vertical direction of the picture and for the values 
t = 0, i, - the same set of lightcones in the planes spanned by r and cf> is plotted. 
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(b) t = 1/a 




(c) t = 



Figure 9. Inside the torus T lightcones tip over as t increases. 
T is represented by the blue stripes and the red curve is given by 
p = 0. The latter becomes a closed null curve at t = - and is a 

' a 

CTC for t > i. This illustration is due to T. Schonfeld 

[2] A. Ori. A class of time-machine solutions with compact vacuum core. Phys. Rev. Lett. 95, 
021101, 2005. 

[3] Geometric analysis of particular compactly constructed time machine spacetimes J. Geom. 
Phys., in press, 2011 

[4] B. O'Neill. Semi-Riemannian Geometry. With Applications to General Relativity. Academic 
Press, Inc., 1983. 

[5] A. Ori. Formation of closed timelike curves in a composite vacuum-dust asymptotically flat 
spacetime. Phys. Rev. D 76, 0440021, 2007. 

[6] S. W. Hawking and G.F.R. Ellis. The large scale structure of space-time. Cambridge Univer- 
sity Press, 1973. 

[7] A. Ori and Y. Soen. Causality violation and the weak energy condition. Phys. Rev. D 49, 
3990, 1994. 

[8] A. Ori and Y. Soen. Improved time-machine model. Phys. Rev. D 54, 4858, 1996. 
[9] T. Schonfeld. On the visualization of geometric properties of particular spacetimes. Master's 
thesis, TU Berlin, 2009. 

23 



[10] A. S. Petrov. Einstein Spaces. Pergamon Press, 1969. 



2 Department of Mathematics, Technische Universitat Berlin, Str. d. 17. Juni 136, 
10623 Berlin, Germany 

E-mail address: scherfneramath.tu-berlin.de 



24 



